Computer control of gene expression: Robust 
setpoint tracking of protein mean and variance 

using integral feedback 
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Abstract — Protein mean and variance levels in a simple 
stochastic gene expression circuit are controlled using propor- 
tional integral feedback. It is shown that the protein mean level 
can be globally and robustly tracked to any desired value using 
a simple PI controller that satisfies explicit sufficient conditions. 
Controlling both the mean and variance on the other hand 
requires the use of an additional control input, chosen here as 
the mRNA degradation rate. Local robust tracking of mean and 
variance is proved to be achievable using multivariable PI control, 
provided that the reference point satisfies necessary conditions 
imposed by the system. Even more importantly, it is shown that 
there exist PI controllers that locally, robustly and simultaneously 
stabilize all the equilibrium points inside the admissible region. 
Simulation examples illustrate the results. 

Index Terms — Gene expression network; moment control; PI 
control; absolute stability. 

I. Introduction 

Synthetic biology is an emergent field of 
biology/biotechnology in which living cells are genetically 
modified in order to achieve new functions. Biological 
circuits of interacting genes and proteins have been 
successfully introduced into living cells to implement various 
functioning modules such as oscillators Q], |2j, switches (3), 
logic gates Q among others. Inspired by electronic circuit 
design, synthetic biologists aim to design functioning modules 
that can be put together in various configurations to build 
biological devices with designed function. In spite of many 
successes, circuits with more than a few interacting genes can 
rarely built reliably and predictably. For this reason, when it 
comes to regulation, building effective biosynthetic control 
circuits remains a big challenge. 

Recently it has been shown that more complex regulation 
becomes possible by moving control functions outside the cells 
and relegating these to a digital computer [5|. This 'in-silico' 
regulation of living cells also allows the implementation of 
controllers with more accuracy and higher speeds than what 
is possible within the environment of the cell. Using synthetic 
biology methods, yeast cells were genetically engineered so 
that their gene expression was responsive to light signals. At 
the same time, flow cytometry and microscopy methods were 
used to measure and quantify the resulting protein expression 
levels in real time. These measurements were then fed to 
a computer control system that used a Kalman filter/MPC 
control configuration to generate a control signal consisting 
of a train of light pulses that drove gene expression. Using 
this setup, feedback control of gene expression in living cells 



was successfully demonstrated experimentally j5]. 

Preliminary results on moment control of reaction networks 
using PI control laws have been obtained in (6) where it is 
shown that such controllers can achieve the desired objectives 
for certain simple reaction networks. The problem addressed 
in this paper is slightly different and directly inspired from 
the relevant gene expression network considered in [5 |. In the 
current work, we use a stochastic model of gene expression 
to explore the feasibility of using simple controllers and 
continuous (non-pulsed) control inputs to achieve effective 
genetic control of protein mean and variance. The problems 
addressed in the paper are beyond those of [6] since existence 
and characterization of positive PI controllers, accounting for 
the presence of nonlinearities, are discussed. It is shown that 
local regulation of protein mean levels to any desired value 
is achievable by acting on the DNA transcription rate using 
proportional-integral feedback. The PI controller is also shown 
to be locally robust and exact regions in the controller param- 
eter space that maintain this local robustness are derived. The 
positivity requirement for control, ignored in |6|, introduces 
a static nonlinearity in the feedback system and must be 
considered to rigorously characterize global stability of the 
controlled system. Using absolute stability theory [7], |8| and 
the Popov criterion [9], it is proved that global asymptotic 
stability is achievable when the controller gains satisfy very 
mild sufficient conditions. 

Next, we show that using a second control input that 
controls mRNA degradation, a multivariable PI feedback con- 
troller can be designed so that any desired protein mean 
and variance setpoints that lie within a certain admissibility 
region defined by the system can be locally and robustly 
tracked. More importantly, it is also shown that there exist 
common multivariable PI controllers that locally and robustly 
stabilizes all the equilibrium points in the admissible region. 
Numerical simulations finally demonstrate the effectiveness of 
the designed genetic control systems. 

Outline: In Section [TTJ the general framework and the main 
problem, i.e. the control of the moments related to the master 



equation, of the paper are introduced. Sections III and IV 



respectively address the problems of controlling the mean 
and the variance of the number of proteins in a simple gene 
expression network. Examples are treated in Section [V] 

Notations: The notation is standard. Given a random vari- 
able X, its expectation is denoted by E[X]. For a square matrix 
M, Sym[M] stands for the sum M + M T . Given a vector 



v € K n , the notation diag(u) stands for a diagonal matrix 
having the elements of v as diagonal entries. 

II. Problem Statement 

A. General framework 

Let us start with the general stochastic formulation where 
N molecular species Si, ... , Sn interact with each others 
through M reaction channels R\,.. . ,Rm- Assuming homo- 
geneous mixing and thermal equilibrium, the time evolution 
of the random variables Xi(t), . . . , Xjy(t) associated with the 
population of each species can be described by the so-called 
Chemical Master Equation (CME), or Forward Kolmogorov 
equation, given by 



M 



^ [w k {x - s k )P{x - s k , t) - Wk(x)P(x, t)] 
k=l 



where s k is the stoichiometry vector associated with reaction 
Rf. and Wk the propensity function capturing the rate of the 
reaction Rk- The variable x is the state-variable and P(x,t) 
denotes the probability to be in state x at time t. 

Based on the CME, dynamical expressions for the first- and 
second-order moments may be easily derived and are given by 

"§3 = SEK*)], 
dE[XX T 



dt 



= SE[w(X)X T ] + E[w(X)X T ] T S T 
+Sdmg{E[w(X)]}S T 



(1) 



3 NxM 



is the stoichiometry 



w li] 



the propensity 



where S := Isi 
matrix and w(X) := \i 
vector. 

According to the structure of the propensity functions 
w(X), the above set of equations may suffer from well- 
posedness and closedness problems. This is explained by 
the fact that computing a moment of a certain order may 
require the knowledge of higher order moments. It is therefore 
not possible, in this case, to describe the evolution of the 
the first moments by a finite set of ordinary differential 
equations. Some approximation schemes have been proposed 
in the literature to overcome this problem, see e.g. |10|-|12|, 
where it is proposed to approximate higher order moments as 
functions of lower order ones. 

B. Ajfine propensity case and gene expression 

In the affine propensity case, i.e. w(X) = WX + wq, 
W G R MxN , w G M M , things turn out to be much 
more convenient. In this case, the moments equations can be 
reformulated as 



dE[X] 
dt 



SWE[X] + Sw , 

Sym[SWS] + Sdiag(WE[X] + w )S T 



(2) 

where S := E[(X - E[X])(X - E[X}) T ] is the covariance 
matrix. It is immediate to see that when the matrix SW 



is Hurwitz, the mean and variance trajectories exponentially 
converge to the equilibrium points X and S given by 

X = -(SW^Swq 



£ = 



.W 1 S 1 s 



Sdi&g(WX + w )S T e SWs ds 



(3) 
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R\ : <j> mRNA 

R 2 : mRNA ^> cf> 

fc 

Rz : mRNA protein+mRNA 
i?4 : protein ^> <j> 



Fig. 1. Simple gene expression network. 

Considering now the simple gene expression circuit depicted 
in Fig. [TJ which involves two species, i.e. mRNA (Si) and 
protein (S2) species, and four reaction channels R\,. .. ,i?4. 
Let Xi(t) and X2(t) be random variables describing the 
number of mRNA and protein molecules respectively. The 
stoichiometry matrix associated with the gene expression 
network is given by 



S = 



-1 
1 







and the (affine) propensity vector by 

w(X) = [k r J r X! k v X x JpX 2 \ 
In vector form, equations Q rewrite 



(4) 



(5) 



xi(t) 
x 2 (t) 



± 3 (t) 
±i{t) 
x 5 {t) 













Xl(t) 

x 2 (t) 



x 3 (t) 

Xi{t) 

x 6 (t) 



B, 



B„ 



k T (6) 



where the state variables are defined by 

= E[X] and 
and the system matrices by 



X 3 Xi 

X4 X5 



A,. 



~7r 

-2 lr 
k p 






-7p. 



A n 





-{lr +7p) 



fop 






(7) 

where k r > is the transcription rate of DNA into mRNA, 
j r > is the degradation rate of mRNA, k. p > is the 



translation rate of mRNA into protein and j p > is the 
degradation rate of the protein. 

Property 1: System (|6|-((7]l is asymptotically stable. 

Property 2: The equilibrium point of system <|6]>-([7j is 
given by 



7r 



kp kf 



h. )■ 
It 



kpkrjjp + k p +7 r ) 

7r(7p+7r)' ^ 7p7r(7p+7r) 

III. Mean control 



(8) 



The objective of the current section is to give a clear 
picture of the mean control of the number of proteins using 
a simple positive PI controller, i.e. a PI controller generating 
nonnegative control inputs. The considered control input is 
the transcription rate k r which has been shown to be possibly 
externally actuated using, for instance, light-induced transcrip- 
tion [5|. It is shown in this section that a positive PI control 
law allows to achieve global and robust output tracking of the 
mean number of proteins. 

A. Preliminaries 

The considered system consists of the interconnection of 
the following restriction of system d6|: 



ii(t)' 

±2{t) 



= A P 



xi(t) 
x 2 {t) 



B e u(t) 



(9) 



and the positive PI control law 

u(t) = if - x 2 (t)) + k 2 I [//* — X2(s)]ds ) (10) 



where ^t* is the mean number of protein to track and the 
scalars k\ , k 2 the gains of the controller. The on-off nonlin- 
earity <p(u) := max{0, u}, see e.g. fl3| , is considered in order 
to impose the positivity of the control input, i.e. the positivity 
of k r = u. 

Property 3: Given a constant reference ^t* > 0, the equi- 
librium point of the system |9|-([To| is given by x* 2 = /Lt* and 

* _ A**7 P * _ A^7 P 7r „ _ u* 
1 ~ ~k~' U ~ ~~k ' ~~ IT 

n>p njp n>2 



(11) 



where 7* is the equilibrium value of the integral term. 

As expected, the presence of the integrator allows to rule out 

any steady state error in the constant reference case. 

B. Local stabilizability, stabilization and output tracking 

Since the equilibrium control input u* and the reference p,* 
are simultaneously positive, the nonlinearity (p(-) is not active 
in a sufficiently small neighborhood of the equilibrium point 
(jTTJi. The on-off nonlinearity can hence be locally ignored 
and the local analysis performed on the corresponding linear 
system. Assuming first that the system parameters k p ,j p ,j r 
are exactly known, the following result on local nominal 
stabilizability and stabilization can be obtained: 

Lemma 4: Given system parameters &p,7p,7r > 0, the 
system |9]) is locally stabilizable using the control law (lOi. 



Moreover, the equilibrium point ( 1 1 1 of the closed-loop system 
(|9|-([T0| is locally asymptotically (exponentially) stable if and 
only if the conditions 



ki > 



7p7r 



and k 2 > 



(12) 



hold. A 
Proof: The local augmented system gathering the linear 
dynamics of the system ([9]) and the controller ( 10 1 is given by 

-k% k 2 





"ii(t)" 




x 2 (t) 




i(t) 





-It 







-7 P 
-1 





~Xl(t)~ 




V 




x 2 (t) 


+ 







J(t) 




1 



[i* (13) 



where I is the integrator state of the controller. Local stabiliz- 
ability is then equivalent to the existence of a pair (fci, k 2 ) € 



M 2 such that the state matrix of the augmented system dl3j 
is Hurwitz, i.e. has poles in the open left-half plane. The 



Routh-Hurwitz criterion yields the conditions (|T2ji that define 
a nonempty subset of the plane (ki,k 2 ). System (|9| is hence 
locally stabilizable using the PI control law ( flC)] ) for any triplet 
of parameter values (fe p ,7 r ,7 P ) € K>o- As a consequence, 
the closed-loop system is locally asymptotically stable when 
the control parameters are located inside the stability region 
defined by the conditions 



12) 



In order to extend the above result to the uncertain case, we 
assume here that the system parameters (fc p ,7 p ,7 r ) belong to 
the set 



(0,fc+] x [ lp ,oo) x [ 7r ,oo) 



(14) 



where the parameter bounds k p , 7~ and 7^ are real positive 
numbers. We then obtain the following result: 

Lemma 5: The system |9]l with uncertain constant param- 
eters (fe p ,7j,,7r) S Vn is robustly locally asymptotically (ex- 



ponentially) stabilizable using the control law ( 10 1. Moreover, 
the equilibrium point (Hi of the closed-loop system |9]h 10 1 
is locally robustly asymptotically (exponentially) stable if and 
only if the conditions 



fci > 



IT 



7p +7r 



and k 2 > (15) 



hold. 



Proof: Define the lower bound function for k\ by 

f(x,y,z) := — ^— - — , x,y,z,k 2 > and let 
y + z x 

f := sup f(x,y,z). Simple calculations show that 

f(x, y, z) is increasing in x and decreasing in y, z over 
(x,y,z)_ e TV Hence, we have / = f{k+,j~,j~) and 
ki > / implies that fci > f(x,y,z) for all (x,y,z) e V^. 
This concludes the proof. ■ 

C. Global stabilizability, stabilization and output tracking 

Local properties obtained in the previous section are gen- 
eralized here to global ones. Noting first that the nonlinear 
function ip(-) is time-invariant and belongs to the sector [0, 1], 
i.e. < ip(x)/x < 1, x € ffi, stability can then be analyzed 
using absolute stability theory fl4[ and an extension of the 



Popov criterion (9), fl4) for marginally stable systems fl5) , 
p6) . We have the following result: 

Theorem 6: Given system parameters fcp,7 P ,7r > 0, then 
the equilibrium point (Hi of the closed-loop system |9]i-(10l 
is globally asymptotically stable if the following conditions 

fc 2 

(16) 



k 2 > and fci > 



1p 



hold. A 
Proof: Following the absolute stability paradigm, the 
closed-loop system (|9|)-([T0|) is rewritten as the interconnection 
of the marginally stable LTI system 

Vk 2 ) 



H{s) 



k p {kis 



(17) 



s(s + 7r)(s + 7 P ) 

and the static nonlinearity </?(•). We assume in the following 
that k 2 > 0, which is a necessary condition for local asymp- 



totic stability of the equilibrium point (111. Since lim sH(s 



p«<2 



> 0, then the Popov criterion |9|, 



7r7p 

It states that the system Q is abso 
controller 



16 



can be applied. 



utely stabilizable using 



10 1 if there exist (fci,fc 2 ) € R and q > such 



that the condition 



F{ju,q) :=%l[{l + qju)H{ju>)] > -1 



(18) 



holds for all uj 6 K.. In order to check this condition, first 
rewrite F(joj,q) as 



F(jUJ,q) 



N (uj) , ^(w) 



D(w) D(w) 



(19) 



where 



£>(w) 



= (or 



fcl(7r7p - - k 2(lr + 7p)] 

fciw 2 (7 r + 7 P ) + k 2 (jpj r ~ w 2 )] 
+ 7 2 )(- 2 ' 



7 2 ) 



Since D(w) > for all u € 
equivalent to 



then the condition ( 1 8 1 is 



N (uj) + qN 1 (uj) + D(uj) > 
for all uj € K. Letting := cj 2 , we get 

Z(u>) := (I) 2 + Zi(<ji)w + z (q) > 
for all Q € [0, oo) where 



(21) 



zo(«) 



7 2 7 2 + fc p [7r(7 P fc i - fc 2 + 7pfe(?) ^ 7p^2 



z i(l) = 7p + 7,- + fcp [(fci7p - k 2)q - h + fci7 r 

The problem therefore essentially becomes a positivity anal- 
ysis of the polynomial Z(ui) over [0,oo). By virtue of 
Descartes's rule of signs p7| , if there exists q > such that 
zo(g) > and > 0, then Z does not have positive 

real zeros, and hence the Popov condition is verified. Since 
&2 > 0, it is immediate to see that by choosing a sufficiently 
large q > 0, the coefficient Zo(q) can be made positive. When 
additionally ki^ p — k 2 > holds, then Z\{q) can also be made 
positive in the same way. Therefore, when these conditions on 
the controller gains hold, Zo(q) and Z\{q) both admit positive 



values provided that q > is chosen sufficiently large, proving 
then that the equilibrium point (JTTJ is globally stable. 

To prove global asymptotic stability, it is enough to note 
that since u* > 0, we have <p(u*) — u* and the control input 
equilibrium value does not lie in the kernel of (p. According to 
p6[ , this allows to conclude on the global asymptotic stability 
of the equilibrium point 



Hi. The proof is complete. 



It is immediate to obtain the following extension to the 
uncertain case: 

Lemma 7: Given system parameters (k p ,j p ,j r ) £ V^, then 
the equilibrium point (Hi of the closed-loop system (|9|- 
( |T0] > is globally robustly asymptotically stable if the following 
conditions 



> and ki > 



1p 



(23) 



hold. A 
It is interesting to note that the above results are also valid 
for any static nonlinearity in the sector [0, 1], including for 
instance saturations. It can also be easily adapted to consider 
nonlinear functions inside more general sectors, albeit the 
resulting conditions may be different. 

Note that complementary sufficient conditions for nominal 
and robust stability could have been extracted from the polyno- 
mial Z in (22 1, either by considering the cases fci7 P — fc 2 = 
and fci7 P — fc 2 < 0; or by using Sturm series JTSJ that could 
provide necessary and sufficient conditions for the polynomial 
Z to be positive over uj e [0,oo). These conditions are 
however quite intricate and, for the sake of simplicity, have not 
been retained in the current work. It is indeed quite to difficult 
to draw interesting conclusions from them, unlike conditions 



(20) ( 16 1 which can interpreted as a restriction of the local stability 



conditions ( |12) , 

D. Disturbance rejection 

It seems important to discuss disturbance rejection proper- 
ties of the closed-loop system. Due to the nonlinear term, the 
rejection of constant input disturbances is only possible when 
they remain within certain bounds. This is formalized below: 
(22) Lemma 8: Given system parameters k p ,j p ,j r > 0, the 



control law ( 10 1 globally rejects constant input disturbances 
5 U that satisfy 

(24) 



o u < — — fit 



provided that the controller gains satisfy conditions (16 1. A 
Proof: In presence of constant input disturbances, the 
equilibrium value of the control input is given by 



7p7r 



(25) 



This value needs to be nonnegative in order to be driven by 
the on-off nonlinearity ip, which is the case if and only if 
condition |24) holds. ■ 
The above result readily extends to the uncertain case: 



Lemma 9: Assume (fc p ,7 p ,7 r ) £ Vu, the control law 
(jTOj satisfying conditions ( |23j ) globally and robustly rejects 
constant input disturbances if and only if the condition 



5„ < 



lp lr 



(26) 



is fulfilled. A 
Proof: The proof follows from a simple extremum argu- 
ment similar to the one used in the proof of Lemma [5] ■ 
It seems important to point out that the sets of admissible 



perturbations defined by ( 24 1 or ( 26 1 do not depend on the 
choice for the controller gains. They do however depend on 
the mean reference value jit*, which is expected since small 
/i*'s yield small control inputs that are more likely to be 
overwhelmed by disturbances. 

E. Concluding remarks 

The coefficient of variation C v := a* /(J,*, defined as the 
ratio of the equilibrium values for the standard deviation and 
the mean number of proteins, is given in the current setup by 



C u = 




(27) 



This shows that the equilibrium standard deviation depends 
on the desired mean value, and thus we have no control over 
it. Since the variance automatically increases as the mean 
increases, this motivates the aim of controlling the variance 
in order to keep it a reasonably low level. 

IV. Mean and Variance control 

As discussed in the previous section, acting on k r is not 
sufficient for controlling both the mean and variance equilib- 
rium values. It is shown in this section that variance control 
can be achieved by adding the second control input j r = U2- 
Fundamental limitations of the control system are discussed 
first, then local stabilizability is addressed. 

A. Fundamental limitations 

Let us consider in this section the control inputs k r = ui 
and j r = U2- It is shown below that there is a fundamental 
limitation on the references values for the mean and variance. 

Proposition 10: The set of admissible reference values 
(jU*, al) is given by the open and nonempty set 

A:=^(x,y)eR 2 >Q : x<y< (l + ^jx^ (28) 

where k p , -f p > 0. A 
Proof: The lower bound is imposed by the coefficient of 
variation which gives 



oi= l 



lr + lp 



(29) 



The upper bound is imposed by the positivity of the unique 
equilibrium control inputs values given by 



1p 



l L * u 2i 



'lp 



(30) 



which are well-posed since a 2 — /i* > according to the coef- 
ficient of variation constraint. The second equilibrium control 

input value u* 2 is positive if and only if erf < I 1 



7> 
lp 



which in turn implies that u\ is nonnegative as well. The proof 
is complete. ■ 
The lower bound obtained above remains valid when k p or 
7 P are chosen as second control inputs. The factor of the 
upper-bound however changes to l + fc p /7 r when 7 p = u 2 , or 
becomes unconstrained when k p = U2- Note however that the 
upper bound on the variance is not a strong limitation in itself 
because we are mostly interested in achieving low variance. 

Note also that since the lower bound on the achievable 
variance is independent of the controller structure, it is hence 
pointless to look for advanced control techniques in view of 
improving this limit. A positive fact, however, is that the 
lower bound is fixed and does not depend on the knowledge 
of the parameters of the system. This potentially makes low 
equilibrium variance robustly achievable. 



B. Problem formulation 

Considering the control inputs k r = u\ and j r 
system (|6]l can be rewritten as the bilinear system 



i«2, the 



h 
i 2 



—U2X1 + ui 

k p xi - j p x 2 
u 2 xi — 2u 2 x 3 
k p x 3 - j p x 4 - 
k p xi + j p x 2 + 



t- Ui 
U2X4 

^k<pX 4 



27p£5 



(31) 



X2 
£5 



where 1\ and I 2 are the states of the integrators. The control 
inputs are defined as the outputs of a multivariable positive PI 
controller 

ui = ip {kid + k 2 h + k 3 e 2 + k^l 2 ) 
u 2 = if (fc 5 ei + k 6 Ii + k 7 e 2 + k 8 I 2 ) 



(32) 



where ex :— /i* — x 2 and e 2 := 07 — x 5 . 

Property 11: Assume that fefcs — k 4 k^ 7^ 0, then the 
equilibrium point of the system pT|H 32 1 is unique and given 

by 



* 1p 
x x = —n 



x 2 = l« 
_ lp 



/ i * u 2; u 2 



'lp 



lp + -> 



and 



" 1*1 ' 




' k 2 


k 4 


-1 




. % _ 




k 6 


k& 




u 2 



(33) 
(34) 



Associated with the set of admissible references A, we 
define the set of equilibrium points as 

X* :={(**,/*) SM 7 : (y.,o*)€A}. (35) 

C. Local stabilizability and stabilization 

Since the equilibrium control inputs are positive, the non- 
linearities are not active in a neighborhood of the equilibrium 
point ([33|-([34|. Local analysis can hence be performed using 
standard linearization techniques. The corresponding Jacobian 
system is given by 

ii = A}xi (36) 
where A\ is given in (37 1 with S := /x* — a\. 



The following result states conditions for the Jacobian 
system to be locally representative of the behavior of the 
original nonlinear system: 

Lemma 12: The Jacobian system fully characterizes the 
local behavior of the controlled nonlinear system (31l-(32l 
if and only if the condition k 2 k$ — k^k§ ^ holds. A 
Proof: For the Jacobian system to represent the local be- 
havior, it is necessary and sufficient that A\ has no eigenvalue 
at 0. A quick check at the determinant value 

det(Al) = ^ p k p (k 2 k 8 - k 4 k 6 )(^(k p + j p ) - j p a%) 

yields that the condition k 2 k% — k 4 ka ^ is necessary and 
sufficient for the local representativity of the nonlinear system. 
Note that since (//*, of) £ A, the term /i*(fc p + 7 P ) — 7 p of is 
always different from 0. The proof is complete. ■ 
The local system being linear, the Routh-Hurwitz criterion 
could have indeed be applied as in the mean control case, but 
would have led to very complex algebraic inequalities, difficult 
to analyze in the general case, even for simple controller struc- 
tures. The Popov-Belevitch-Hautus (PBH) stabilizability test 
would not have helped either to conclude on anything useful 
since it does not take into account the controller structure. 
Despite the 'large size' of the matrix A\, it is fortunately still 
possible to provide a stabilizability result using the fact that A* t 
is marginally stabl^jwhen the control parameters fej are set to 
0. This is obtained using perturbation theory of nonsymmetric 
matrices J2Q] . 

Lemma 13: Given any k p ,j p > 0, the bilinear system (31 
is locally asymptotically stabilizable around any equilibrium 
point (33l-(34i using the control law (32i. A 
Proof: The perturbation argument relies on checking 
whether the eigenvalues on the imaginary axis can be shifted 
by slightly perturbing the controller coefficients around the '0- 
controller', i.e. by letting fcj = edi, where e > is the small 
perturbation parameter and cfj is the perturbation direction for 
controller parameter fej. We assume here that both integrators 
are involved in the controller, that is jc^l + |de| > and 
| (£4 1 + \dg\ > 0. To prove the result, let us first rewrite the 
matrix A| as 

8 

A* e = Ao + e^djAj. (38) 



A* t has eigenvalues in the closed-left half plane and those on the imaginary 
axis are semisimple 1191. 



The matrix A n is a marginally stable matrix with a semisimple 
eigenvalue of multiplicity two at zero. Paradoxically, these 
eigenvalues introduced by the PI controller are the only critical 
ones that must be stabilized, i.e. shifted to the open left-half 
plane. From perturbation theory of general matrices [20], it is 
known that semisimple eigenvalues bifurcate into (distinct or 



not) eigenvalues according to the expression |20| 



Xi(e, d) = e£i(d) + o(e), i=l,2 



(39) 



where £i(d) is the 

8 

M(d) := ^2 d * M * with 



( = 1 



eigenvalue of the matrix 



(40) 



Above, v\ , v\ and v\ , are the normalize^ left- and right- 
eigenvectors associated with the semisimple zero eigenvalue. 
It turns out that all Mi's with odd index are zero, indicating 
that the proportional gains have a locally negligible stabilizing 
effect. This hence reduces the size of the problem to 4 
parameters, i.e. those related to integral terms. We make now 
the additional restriction that cLt = d$ = reducing the 
controller structure to one integrator per control channel. The 
matrix M (d) then becomes 



M(d) = tp 



where ip := 



k p d 2 
k p (7^d 2 

- 



—daft* 



-2ai 



* . The semisimple eigen- 

7p(M* - cri) + fcpA** 
values then move to the open left-half plane if there exist 

perturbation directions d 2 ,ds € K, d 2 dg ^ 0, such that M(d) 

is Hurwitz. We can now invoke the Routh-Hurwitz criterion 

on M{d) and we get the conditions 



d 2 d s ip 
2 



7 P M* 



> 



< 




0. 



7 P V (da (7p(m* - °"*) 2 - 2fc p ^*cr^) + d 2 k p ^<jl) 

Since the term ip is negative for all (/i*,cr^) e A, the first 
inequality holds true if and only if d 2 d$ < 0, i.e. perturbation 
directions have different signs. The second inequality can be 
rewritten as 



d 2 > da 



(41) 



Choosing then d s < 0, there always exists d 2 > such 
that the above inequality is satisfied, making thus the matrix 
M(d) Hurwitz. We have hence proved that for any given pair 



£ A, there exists a control law (32i that makes the 



corresponding equilibrium locally asymptotically stable. The 
proof is complete. ■ 



2 Normalized eigenvectors verify the conditions v\v\ 
}v~ = vful = 0. 
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(37) 



It is possible to go beyond this result and show that there 
exist semi-global PI controllers: 

Lemma 14: Given any fc p ,7 P > 0, there exists a common 
control law ([32]) that simultaneously locally asymptotically 



stabilizes system ( 3 1 1 around all the equilibrium points in X . 

A 

Proof: To show that there exists a common controller that 
simultaneously makes all the equilibrium points in X* locally 
asymptotically stable, it is enough to prove that there exists a 
pair (d,2,ds) G K 2 , d 2 ds such that the inequality (41 1 is 
satisfied for all (/x*,af) € A. This is equivalent to finding a 
finite e?2 > satisfying 



d 2 > d 8 2 



sup 



7p(/«* - at) 



2\2 ' 



kpfji^o"^ 



(42) 



Standard analysis allows to prove that 



f 7 P (/x, -al) 2 \ _ 



SUp s 

„<rl)eA I «pM*cr; 



7p 



G (0, 1) (43) 



which shows that by simply choosing the directions d$ < 
and d 2 > 0, the matrix M(d) becomes then Hurwitz for 
all (/i*.cr 2 ) e A This therefore implies the existence of 
a common control law ( |32] > that locally and asymptotically 
simultaneously stabilizes all the equilibrium points in X*. ■ 

V. Examples 

For simulation purposes, we consider the normalized ver- 
sion of system pTj ) similarly as in |5|: 

-(7° + "2)^l(^)^ 

j p (xi(t) - x 2 {t)) 



x x {t) 
x 2 (t) 
x 3 (t) 
I 4 (t) 

x 5 (t) 



•Ui(t) 

- 2(7° 



(_7r + 7p)S3(*) - (7r + "2 + 7 p)z 4 (£) 

[xi{t) + x 2 {t) + 2(a - l)x 4 (t)} - 2 lp x 5 (t) 



7p 



where %(t) = 7 ° + bui(t), 7 ° = 0.03, 7 P = 0.0066, 
b = 0.9587, k p = 0.06 and a = 1 + fcp/( 7 ° + 7p ). The system 
has been normalized according to basal levels for transcription 
rate k® and degradation rate 7 J?. In the absence of control 
inputs, i.e. ui = and U2 = 0, the system converges to 
the normalized equilibrium values x* = 1, i = 1, . . . , 5. The 
parameter values are borrowed from |5]. 



A. Mean control 

The considered PI controller parameters computed using 
loop shaping are k\ = 0.01 and k 2 = 0.0007. Simulations 
yield the trajectories of Fig. [2] and Fig. [3] We can see that, 
as expected, the proposed controller achieves output tracking 
for different references and in presence of constant input 
disturbances. 
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Fig. 2. Trajectories of the mean number of proteins for different reference 
values. 



B. Mean and variance control 

In this case, a PI controller with gain k\ = 1, k 2 = 0.007, 
k 3 — —0.2 and k 4 = —0.0014 is considered. The normalized 
achievable minimal variance is given by 



X + 7p 

7r + 7p + K 



-11. 



(44) 



The response of the controlled variance according to changes 
in the reference value is depicted in Fig. [4] where we can 
see that the variance tracks the desired value quite well. In 
order to avoid oscillations, the changes in the reference values 
follow a ramp. It seems also important to point out that when 
the reference point changes, due to the coupling between the 
mean and variance, the mean value changes as well, but this 
is immediately corrected by the mean controller. 



5 



Finally, the implementation on the real process is under way. 
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Fig. 3. Trajectories of the mean number of proteins in response to constant 
input disturbance. 
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Fig. 4. Response of the controlled variance according to changes in the 
reference. 



VI. Conclusion 

Controlling the mean and variance of the number of proteins 
in a simple gene expression circuit using PI controllers has 
been shown to be achievable. Interestingly, PI controllers have 
been proved to be sufficient for respecting most of the inner 
(system) and outer (design) constraints, and to be global in 
the sense that a single PI can locally stabilize all the possible 
equilibrium points. The very same results have been extended 
to uncertain systems. 

Future works will be devoted to a better characterization 
of the stability domain of the controlled variance dynamics 
(invariant set, global stability, etc), the derivation of advanced 
PI controllers to improve transient behavior and addressing 
implementation issues. The generalization of this idea to the 
non-affine propensity case and the control of higher order 
moments are also important problems that will be considered. 
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